Modified f(G) gravity models with curvature-matter coupling 
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A modified f(G) gravity model with coupling between matter and geometry is proposed, which 
is described by the product of the Lagrange density of the matter and an arbitrary function of 
the Gauss-Bonnet term. The field equations and the equations of motion corresponding to this 
model show the non-conservation of the energy-momentum tensor, the presence of an extra-force 
acting on test particles and the non-geodesic motion. Moreover, the energy conditions and the 
stability criterion at de Sitter point in the modified f(G) gravity models with curvature-matter 
coupling are derived, which can degenerate to the well-known energy conditions in general relativity. 
Furthermore, in order to get some insight on the meaning of these energy conditions, we apply them 
\j—, • to the specific models of f(G) gravity and the corresponding constraints on the models are given. 

In addition, the conditions and the candidate for late-time cosmic accelerated expansion in the 
modified f(G) gravity are studied by means of conditions of power-law expansion and the equation 
of state of matter less than 
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I. INTRODUCTION 



According to recent observational data sets[l|, |2J, our current universe is flat and undergoing a phase of the 

ly-, ' accelerated expansion which started about five billion years ago. In principle, this phenomenon can be explained by 

either dark energy (see, for instance, Ref.[3| for reviews), in which the reason of this phenomenon is due to an exotic 

cn ! . . . n 

component with large negative pressure, or modified theories of gravity 4] . Unfortunately, up to now a satisfactory 
answer to the question that what dark energy is and where it came from has not yet to be obtained. Alternative 
to dark energy, modified theories of gravity is extremely attractive, such as f(R) gravity (see, for instance, Ref. 5] 



for reviews), here f(R) is an arbitrary function of the Ricci scalar R. Cosmic acceleration can be explained by f(R) 
gravity |6|, and the conditions of viable cosmological models have been derived in [7|. 



A general model of f(R) gravity has been proposed in Ref.[8|, which contains a non-minimal coupling between 
geometry and matter. This coupling term can be considered as a gravitational source to explain the current acceler- 
ation of the universe. As a result of the coupling the motion of the massive particles is non-geodesic, and an extra 
force, orthogonal to the four- velocity, arises. Different forms for the matter Lagrangian density L m , and the resulting 
extra-force, were considered in [9( , and it was shown that more natural forms for L m do not imply the vanishing of the 
extra- force. The implications of the non-minimal coupling on the stellar equilibrium were investigated in |10|], where 
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constraints on the coupling were also obtained. An inequality which expresses a necessary and sufficient condition to 
avoid the Dolgov-Kawasaki instability for the model was derived in [ill ]. However, a more general model, in which the 
coupling style is arbitrary and the Lagrangian density of matter only appears in coupling term, has been proposed 



in Ref. 12], i.e., the so-called the generalized f(R) gravity with arbitrary coupling between matter and geometry. In 
this class of models the energy-momentum tensor of the matter is generally not conserved, and the matter-geometry 
coupling can induce a supplementary acceleration of the test particles. Moreover, the energy conditions and the 



Dolgov-Kawasaki criterion for the model have been derived in Ref. 13(, which are quite general and can degenerate 
to the well-known energy conditions in GR and f(R) gravity with non-minimal coupling and non-coupling as special 
cases. 

Another interesting alternative modified theory of gravity is the modified Gauss-Bonnet gravity, i.e., f(G) gravity, 
where f(G) is a general function of the Gauss-Bonnet (GB) term[14j, |l5|. At present specific models of f(G) gravity 
have been proposed to account for the late-time cosmic acceleration [15|, Il6| . and the respective constraints on the 



parameters of the models have also analyzed in Ref. lfjj]. For more general forms of f(G), the most crucial condition 

to be satisfied is d 2 f/dG 2 > 0, which is required to ensure the stability of a late-time de Sitter solution as well as the 

I I 
existence of standard radiation/matter dominated epochs 17], and solar system constraints on f(G) gravity models 

have been studied in Ref.[18j. Recently the energy conditions in f(G) gravity have been also discussed [19j, but they are 

only adapted to /(G) gravity without coupling between matter and geometry. Hence, in this paper the /(G) gravity 

models with curvature-matter coupling will be proposed, and some relevant issues, such as the energy conditions, the 

stability criterion and the conditions for late-time cosmic accelerated expansion, will be studied. 

This paper is organized as follows. In Section 2, the /(G) gravity models with curvature-matter coupling are 

proposed, which is here called the modified /(G) gravity. And some fundamental elements of the modified /(G) 

gravity are given. In Section 3, the well-known energy conditions, namely, the strong energy condition (SEC), the 

null energy condition (NEC), the weak energy condition (WEC) and the dominant energy condition (DEC), will be 

derived. Concretely, the two models are applied to the weak energy condition in order to understand the meaning of 

these energy conditions. Furthermore, we will study stability criterion at the de Sitter point, by which the parameters 

in the specific model in the modified f(G) gravity can be constrained in Section 4. In addition, by using the conditions 

of power-law accelerated expansion and the equation of state of matter less than — i, the conditions for late-time 

cosmic accelerated expansion in the modified f(G) gravity are discussed in Section 5. Conclusions and discussions on 

our work are given in the last section. 

II. FIELD EQUATIONS IN THE MODIFIED F(G) GRAVITY 

As we know, the modified Einstein-Hilbcrt action[20| is, 

S = Jd^V-g[^^ + L m ], (1) 

in which k = 8ttGn, Gn is the gravitational constant, R — R(g flu ) is the Ricci scalar, and L m is the Lagrangian 
density of matter. The Gauss-Bonnet invariant is defined as G = R 2 — AR^R^ + R^^aR^^ {R^v and R^^ a are 
the Ricci tensor and the Riemann tensor, respectively). 

Below, we consider f(G) gravity with curvature-matter coupling, which is here called the modified f(G) gravity. The 



Lagrangian density of matter only appears in coupling term and the action is given by 



S 



R 



<TxV=g{-+[l + \f(G)]L m }, 



(2) 



where we have chosen k = 8ttGn = 1, which will be adopted hereafter. Varying the action ^ with respect to the 
metric g^ u yields the field equations: 

Riiv — ^Q^vR = Tfiu + 2XL m (— 2FRR flu + AFR^R U ^ — 2FR fl £ IJ < ; Rl' 7 '' — 'iFR^ c , v R> (T 

+2i?V M V„F - 2Rg^ 2 F - 4R« V,V C F - ARf, V M V ? F (3) 

+4i?^V 2 F + 4g Mt/ i?« CT V e V CT F - 4i? M?t/CT V«V CT F), 

where F = F(G) = df(G)/dG, and T M „ is the energy-momentum tensor of matter, which is defined as: 

2 6(y/=gL m ) 



-*- LLV 



(4) 



-g 8g» v 

By assuming that the Lagrangian density L m of the matter depends only on the metric tensor components, and 
not on its derivatives, we obtain T^ v = L m g^ v — 2dL m /dg liV . By taking the covariant derivative of Eq.Q, with the 
use of the Bianchi identities, V^G^ = ( here G^ v is the Einstein tensor ), we can obtain the following relation: 

VT M „ = 4AL m [V^(iUV - 2R%R vi + R^R^ + 2g^g^R tiav0 R^) 

+F(R IXV V»R + y^RWR - 2R u( y»Rl - 2Rfy*R vi + R^WR^ 

+R^ a ,V^Ri^ + 2g^g^R ia V^R^ p + 2g^g^R mv ^R i<T ) + RV „UF 
-RUV U F - V M V„FV^i? - 2R fi „WnF + 2i?« V^V^V ? F + 2V^V e FV^i?« 
+2i?«DV ? F + 2V M V e .FVi$ - 2g tlu R^V»V ( y a F - 2g^V i V a FV^R^ 
+2g^g^R (iauP V^ i V <T F + 2.g Q «/ ff V ? V CT FV^i? MQ ^], 

from which we see that the conservation of the energy-momentum tensor of matter is violated due to the coupling 
between matter and geometry. 

For the perfect fluid, the energy-momentum tensor is 



T^v = {p + p)u^u v + pg^ v , 

where u^ is the four- velocity, and satisfies the conditions u^u^ = 1 and u^u^ v — 
of Eq.® can be given as 

VT M „ = ( P + p)g^u v V„u» - V„ P (81 - u v u x ). 



(6) 



Thus, the covariant derivative 



(7) 



By imposing the condition of the conservation of the matter current, \7 u (pu u ) = 0, and with the use of the identity 



u"V„u^ = 



d 2 x" 



rtx^Q 



we have the equation of motion of a test particle in the model as 
Du^ du^ „„ „ , d 2 x fJ ' 



(Is 



du* 1 
ds 



T\u»u 



v„,\ 



ds 2 



Kx^ x 



f- 



where the extra-force / M has the following expression 

1 



f 



p + p 



\y»T lu ,g* x + V v p( S r-u v u' 1 )]. 



(8) 



(9) 



By substituting the relation ([5]) into Eq.Q, we can point that due to the presence of the coupling between matter 
and geometry, the motion of the massive particles is non-geodesic, and the extra-force / M is not orthogonal to the 
four-velocity m m , that is, there exists an angle between the extra-force /^ and the four-velocity u^, and f^u^ =/= 0. 



III. ENERGY CONDITIONS IN THE MODIFIED F(G) GRAVITY 
A. The Raychaudhuri Equation 

The energy conditions arise when one refers to the Raychaudhuri equation for the expansion [2l| . Under these 
energy conditions, one allows not only to establish gravity which remains attractive, but also J» keep the demands 
that the energy density is positive and cannot flow faster than light. Below, following Ref.pjj we simply review 
the Raychaudhuri equation which is the physical origin of the null energy condition(NEC) and the strong energy 
condition(SEC) 
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In the case of a congruence of timclikc geodesies defined by the vector field u M , the Raychaudhuri equation is given 
by 

Hf) 1 

-j- = —P - v^ + uvu;^ - R^uW, (10) 

dT 6 

where 11^,6, cr^„ and uj^ v are the Ricci tensor, the expansion parameter, the shear and the rotation associated with 
the congruence, respectively. While in the case of a congruence of null geodesies defined by the vector field k^, the 
Raychaudhuri equation is given by 

HO 1 

-j- = -jr* 2 - ^a^ + uj^lu^ - R^WW. (11) 

dT I 

From above expressions, it is clear that the Raychaudhuri equation is a purely geometric statement and independent 

of the gravity theory. In order to constrain the energy-momentum tensor by the Raychaudhuri equation, one can use 

the Ricci tensor from the field equations of gravity to make a connection. Namely, through the combination of the field 

equations of gravity and the Raychaudhuri equation, one can obtain physical conditions for the energy-momentum 

tensor. Since a 2 = <j liV <j ,il ' > (the shear is a spatial tensor) and u)^ v = (hypersurface orthogonal congruence), 

from Eqs. (|10[) and (TTTj) . the conditions for gravity to remain attractive (M < 0) are 

R^u»u v > SEC, (12) 

R^k v > NEC. (13) 

Thus by means of the relationship (|12j) and Einstein's equation (i? M „ — -^Rg^ = T M „), one obtains 

R^u'u" = (T„„ - -g^uW > 0. (14) 

If one considers a perfect fluid with energy density p and pressure p, 

Tpv = (p + pjUpUv + pg^ v (15) 

the relationship (1141) turns into the well-known SEC in general relativity, i.e., 

p + 3p>0. (16) 

Similarly, by using the relationship (J13I) and Einstein's equation, one has 

T M „fc"fc" > 0. (17) 

Thus, by considering Eq. (|15[) . the familiar NEC in general relativity can be reproduced as: 

P + P>0. (18) 



B. Energy conditions 



The Einstein tensor resulting from the field equation ([3]) can be written 



^flV ^p,U n9\lV^ -* IIV 5 \}-^) 



where 



T$ f = T llv + 2\L m {-2FRR llu +AFR^R vi -2FR llia ,Ri^-AFR^ v R^ 

+2RV fl V u F - 2Rg^V 2 F - 4i?« V„V 4 F - 4i?| V M V ? F (20) 

+4i?^V 2 F + 4.g^i?« CT V c V CT ^ - AR^ va V^V a F). 

Contracting the above equation, we have 

T eff = T _ A\L m (FG + RDF - 2R fiV VV v F). (21) 

Thus, we can write R^ v in terms of an effective stress-energy tensor and its trace, i.e. , 

R^ = T e J f -\g^T*ff. (22) 

By using the relationship (|12l) and Eg. (|2"2"1) . the SEC can be given as: 

T e JJu^u v - -T eff > 0. (23) 

By using the Eqs. (l2"0)) and (J2"T1) . the SEC in Eq. (l2"3")l can be expressed as 

p + 3p + 8XL m FG - 96XL m H 3 i G + A\L m RUF - 8AL m i? M „V A 'V ,y F > 0, (24) 

where H = a(t)/a(t) is the Hubble parameter. The NEC in the modified f(G) gravity can be expressed as: 

T e JJk^k v > 0. (25) 

By the same method as the SEC, the NEC in Eg. (1231) can be changed into 

p + p + AXL m FG - 6AXL m H 3 i G + ^XL m RDF - ^\L m R^V^V»F > 0. (26) 

Furthermore, by means of Eqs. (|2"0"|) and (|2"2")l . the effective energy density and the effective pressure can be derived as 
follows: 

p eff = 2XL m (FG - 2AH 3 f[ G ) + p, (27) 

p e ff = 2XL m FG - 16XL m H 3 i G + ^XL m RDF - ^XL m R^V»V"F + p. (28) 



Note that the above expressions of the SEC (|24l) and the NEC (|26|) are directly derived from Raychaudhuri equation. 
However, equivalent results can obtained by the transformations p — > p ej ' and p — > p Cj ' into p + 3p > and p+p > 0. 
Thus by extending this approach to p — p > and p > 0, the corresponding DEC and the WEC in the modified f(G) 
gravity can be respectively given as: 

p-p- 32XL m H 3 i G - ^XL m RUF + ^XL m R^V»V v F > 0, (29) 
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p + 2XL m FG - 48H 3 f[ G > 0. (30) 

It is worth stressing that when taking f(G) = 0, the energy conditions in gen eral relativity can be reproduced. 
Furthermore, by defining the deceleration, jerk, and snap parameters as|23j 

1 a . 1 a la 

we have 

H = -H 2 (l +q),H = H 3 {j + 3q + 2), H = H 4 (s - 2j - 5g - 3), (32) 

and the GB term is given by 

G = 24H 2 (H 2 +H). (33) 

In addition, we consider L m = —p, we can rewrite Eqs. (|2"T)) and (|2"8)l as follows: 

p e // = p + 48XH 4 p[qf + 24H 4 (2q 2 + 3q + j)f"}, (34) 

p eff = p + 4 8 \ H i p [ q fi + 4H 4 (-2q 3 -8q 2 + q- 2j - 6qj + s + 3)/" + 96H 8 (2q 2 + 3q + j) 2 f"}. (35) 

Hence, the energy conditions, i.e., the SEC, NEC, DEC and WEC can be rewritten as: 

p + 3p + 48XH 4 p[qf + 24H 4 (2q 2 + 3q + j)f"] + lUXH 4 p[qf + 4H 4 {-2q 3 

-8q 2 +q-2j- 6qj + s + 3)f" + 96H s (2q 2 + 3q + iff"] > 0, (36) 

p + p + 48XH 4 p[qf + 24H 4 (2q 2 + 3q + j)f"] + 48XH 4 p[qf + 4H 4 {-2q 3 

-8q 2 +q-2j- 6qj + s + 3)f" + 96H 8 (2q 2 + 3q + j) 2 /'"] > 0, (37) 

P-P + 48XH 4 p[qf + 24H 4 (2q 2 + 3q + j)f"} - 48XH 4 p[qf + 4H 4 (-2q 3 

-8q 2 +q-2j- 6qj + s + 3)/" + 96H s (2q 2 + 3q + j) 2 f'"} > 0, (38) 

p + 48XH 4 p[qf + 24H 4 (2q 2 + 3q + j)f"} > 0. (39) 

It is worth stressing that when taking L m = p, Eqs. (|2"T| and (|28p change into 

p-ff = p - 48XH 4 p[qf + 24H 4 (2q 2 + 3q + 3 )I"], (40) 

p eff ^ p _ 48AiJ 4 p[q/' + 4H 4 (~2q 3 - 8q 2 + q - 2j - Qqj + s + 3)f" + 96H s (2q 2 + 3q + jff'"}. (41) 

And the corresponding energy conditions are as follows: 

p + 3p - 48XH 4 p[qf + 24H 4 (2q 2 +3q + j)f"} - 144XH 4 p[qf + 4H 4 (-2q 3 

-8q 2 + q-2j- 6qj + s + 3)f" + 96H 8 (2q 2 + 3q + j) 2 /"'] > 0, (42) 



p + p- A8XH 4 p[qf + 2AH 4 (2q 2 + 3q + j)f"] - A8XH 4 p[qf + 4H 4 (~2q 3 

-8q 2 + q-2j- 6qj + a + 3)/" + 96H s (2q 2 + 3q + 3? f"\ > 0, (43) 

p-p- 48XH 4 p[qf + 24H 4 (2q 2 + 3q + j)f"] + 48XH 4 p[qf + 4H 4 (-2q 3 

-8q 2 + q-2j- 6qj + a + 3)/" + 96H 8 (2q 2 + 3q + j) 2 /'"] > 0, (44) 

p - A8\H 4 p[qf + 2AH 4 (2q 2 + 3q + j)!"] > 0. (45) 

C. Energy Conditions for specific f(G) models 

In order to exemplify how to use the energy conditions to constrain the f(G) theories of gravity below, we study 
the realistic models of f(G) gravity, which have been found to reproduce the current acceleration [15j,[l6|: 

f 2 (G)=a 3 G n (l + b 3 G m ), (47) 

where 01, 02, 03, 61, 62, &3, n, m are all constants. Since there has been no reliable measurement for the snap parameter 
(s) up to now, we only focus on the WEC (|39[) and (|45|) in this particular case. 

Since the inequalities are too complicated to find exact analytical expressions, so we have to take some specific 
values of the parameters, such as a\ = b\ = — 1, 02 = 2, 62 = b 19]. Also, when G — > ±00 or G — > 0~, the model (jT7|) 
can be changed in the form /(G) ~ aG™[16|, which means 03 = a, b 3 = 0. The above two models can be rewritten as: 

/^) = -g^V (48) 

h{G) = aG n . (49) 

(1) When L m = —p, we can get the corresponding WEC as follows: 

p{l + A8XqH 4 ^t^±^ + U52XH\2q 2 + 3q + j) 

G n - 2 n[2G n (n + l)-(n-l)b}(-2 + b) ^ 
x (2G „ + b)3 >^° Uorfi{G)), (50) 

p[l + A8XqH 4 anG n - 1 + 1152XH 8 {2q 2 + 3q + j)an{n - 1)G"~ 2 ] > {forf 2 {G)). (51) 

After a series of simplification, taking some present values Hq = 70.5[24|, qo = —0.81 ± 0.14 and jo = 2.16^Q75 25], 
and A = 1, we can obtain the restrictions on the parameters n, b and a, which satisfy the WEC in Eqs. (|50| and (|51|) . 
respectively (see Figs.l and 2). Fig.l shows that the WEC p e ^ > can give the constraints on the parameters n 
and b in the fi(G) model, i.e., < b < 1 and —2.3 < n < 2.5, but except \n\ < 0.1 owing to the non-continuity of 
peff > 0. Similarly, Fig. 2 illustrates that only when a > and 7.5 < n < 10, the WEC p e ^ > can be satisfied in 
the / 2 (G) model. 



(2) When L m — p, let us consider a perfect fluid composed of non-relativistic or relativistic particles with constant 



barotropic equation of state (EoS) p = (7 — l)p, < 7 < 2 is a constant relating to the EoS by u> — 7 — 1 26l428| . the 
corresponding WEC changes into: 

[ G n - 1 (-nb + 2n) 



p{l-A8\qcuH 4 



(2G n + by 



ll52XH fi uj(2q z + 3q + j) 



G n - 2 n[2G n (n + 1) - (n - 1)6] (-2 + b) 
(2G™ + 6) 3 



}>0 (forhiG)), 



p[l - A%\qojH i anG n - 1 - 1152XH 8 uj(2q 2 + 3q + j)an(n - l)G n ~ 2 ] > (forf 2 (G)). 



(52) 



(53) 



By the similar discussions to the case of L m = —p, we can obtain the restrictions on the parameters n, b and a, 
which respectively satisfy the WEC in Eqs. (|5"2l and (j5"3"]) . and are illustrated in Figs. 3 and 4 (here taking ui = 0.5). 
From Fig. 3 it is easy to see that the constraints on the parameters n and b for the fi(G) model are —0.2 < b < 1 
and —2.3 < n < 2.5, which are nearly the same as the results of the L m = —p, but in the Fig. 4 —10 < a < 10 and 
—2.7 < n < —1.3 for the /2(G) model are quite different from the results of the L m = —p. 




(a) (b) 

FIG. 1: The constraints of WEC (p e ^ > 0) on the parameters n and b for the /i(G) model in (|48)l with A = 1 and L„ 



= -P- 




(a) (b) 

FIG. 2: The constraints of WEC (p > 0) on the parameters n and a for the /2(G) model in (149)1 with A = 1 and L„ 



= -P- 




(a) 



(b) 



FIG. 3: The constraints of WEC (p e ^ > 0) on the parameters n and a for the fi(G) model in ((48} with A = 1 and L m — p = cop 

(w = 0.5). 




(a) 



(b) 



FIG. 4: The constraints of WEC (p e ^ > 0) on the parameters n and a for the /2(G) model in ([49} with A = 1 and L m — p = cop 
(w = 0.5). 

IV. STABILITY CRITERION AT DE SITTER POINT 



171 , we will focus 



Modified gravity must be stable at the classical and quantum level. There are in principle several kinds of 
instabilities to consider, such as Dolgov-Kawasaki criterion in f(R) gravity [29J- Below, following Ref. 
on the stability criterion at dc Sitter point in the modified f(G) gravity. 
In a flat FLRW background with the metric 



ds 2 



-dt 2 + a{t) 2 dXl. i 



(54) 



where a(t) is the scale factor and dX 2 contains the spacial part of the metric. The 00 component of the field equation 
© gives 



3H 2 = 2XL m Gf, G - A8XL m H 3 l G + Pn 



(55) 



Let us first discuss the stability around the de Sitter point in the modified f(G) gravity by neglecting the contribution 
of pressure- less matter p m and radiation p r . The Hubble parameter, H = Hi (at the de Sitter point), satisfies 

3Hf = 2XL m G 1 f, G (G 1 ) - 48AL mJ ff 1 3 /, G (G 1 ), (56) 



** + J&J' S * + dooir. *?7- ^TT - Tirghrm = 0. (58) 
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where G\ = 24Hf, and the relations H\ = and Gi = are used. Considering a linear perturbation £7?i about the 
de Sitter point, Eq. (|55]) gives 

2\L m {2AHlf fiG 5G 1 - f, G 8Gi) = 8Hi(48\L m H*f jGG - 6#i). (57) 

Substituting the relations 6G{H X ) = 24(4Hf5H x + Hf 8 Hi) and 8G(Hi) = 24Hf(8Hi + 4Hi6H x ) into Eq.dSTJ), we 
obtain 

96 ffj t /, GG -/, G • 1 1 £g_ 

24H?f, GG X [ 192\L m Hlf, GG 24H X 6H 2 f >GG ' 

It follows that the effective mass squared is ( 192AL IF f 2477 6g a? — ), which must be non- negative for 

stability. Therefore, we can obtain 

1 - SXHUgg - 32\L m H 2 f,G > Q (5g) 

192XL m H 1 f t aa 

which is just the stability criterion at the de Sitter point. It follows that if the exact value of Hi and a suitable form 

of L m can be given in the modified f(G) gravity models, then the constraints on the parameters in the specific model 

can be obtained. 

V. THE CONDITIONS FOR LATE-TIME COSMIC ACCELERATED EXPANSION IN THE MODIFIED 

F(G) GRAVITY 

It is known that late-time cosmic accelerated expansion occurs under the conditions of either a power-law 
expansion or the equation of state of matter less than — |. To exemplify how to use these conditions to realize the 
phase of accelerating expansion in the modified f(G) gravity, now we concentrate on the model fi(G) in (|48|) . Thus, 



by means of the action ^ and the energy density L m of perfect fluid[30ll31|. i.e., 

L m = -p=-p a-^ 1+ ^\ (60) 

where uj is the equation of state of perfect fluid and is taken to be a constant, the field equation (|55j) becomes 

3H 2 = -p Q a- 3{1+u) [l - 48\H 2 (H 2 + H)f - 1152H 4 {2H 2 + HH + 4H 2 H)f"\. (61) 

Note that the relation R — 6(2H 2 + H) is used in the derivation of the Eq. ([6"Tl) . 

If assuming the solution of (|0L[) is a = aoi r [13|], then we have H = j, H — — £, j? = |J. Substituting all these 
relations into ([FT]) , we can get the following equation: 

3r 2 = Po «o 3(1+U) ^ 3r(1+W)+2 , 1V , 3 „ x {b 3 (r 1) 

(r - l) 2n r 6 ™ , „ (r - l)"r 3 " „ r(r - l)"r 3 " 
+8- { -!3824™ x ^ / + 13824" x -^ ^ + 5576"nA 

^8n L -fin +An 

(r — ~\ ) n r 3n 

+4 3n+l Tn 2 x _ 576 n nr) l + 2 3n+l x yitf V ) r_ 3 _ ^ + 4 2 A 

J jAn L 

( r _ 1 ^n„3n / _ 1 \n 3n 

+r(3 + nA)] + ^—^ {-3 2 " +1 x 64" Li + 5nA[24" 

-576" l , j ] - 4n 2 A[24" + 2 6 " +1 x 9" [r 'J - 576" l , j ] 

f4n J L +4n jAn ' 

+r[3 2 " +1 x 64"^-^ 24"nA + 576"nA^-^ ]}}. (62) 
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We find six kinds of possible relationships among r, ui and n, namely, r 



_ 2(l-2n) _ 2(1 — in) _ 2(l-6n) 



3(l+w) ' 



"sTT^ 



3(l+u;) ' 



2(1 -8n) 

3(l+w) • 



— ~7i — i — \ — ancl 7* 

3(1+0;) 



2(l-12n) 

3(l+w) ' 



corresponding regions of ui are w < 



3 3" 



Under the condition of a power-law expansion (i.e., r > 1), the 



I | ,• 2(l-2n) 1 8 r 2(l-4n) 1 . r 

" for r = -^t- — t 1 , w < — i — %n for r = -^n — r, to < — 4 — 4n for 



_ 2(l-6ra) 
"" 3(l+w) ' 



< — i — 4^n for r = 



_ 2(l-8n) 



3(l+w) ' 



3(1- 

20, 



3(1- 



< — i — 22n for r = L-, . "^ and w < — i — 8n for r = L-, . T ^ , 

3 3 3(1+oj) 3 3(1+oj) ' 



respectively. Furthermore, by considering the equation of state of matter less than — | (i.e., u> < — |), the relationship 
among r, w and n, condition and candidate for late-time cosmic accelerated expansion are shown in Table 1. The 
candidate for late-time cosmic accelerated expansion can be either the effective quintessence (— 1 < u; < — 4) or the 
effective phantom (u> < —1). 



Relationship 


Condition(r > l,w < — i) 


The effective quintessence 


The effective phantom 


_ 2(l-2n) 
3(l+w) 


n > and n 7^ | 


< n < § 


n>\ 


_ 2(l-4n) 
' 3(l+w) 


n > and n / | 


0<n< i 


n>\ 


_ 2(l-6n) 
' 3(l+u>) 


n > and n/i 


0<n< § 


n>\ 


_ 2(l-8n) 
' 3(l+w) 


n > and n 7^ | 


< n< | 


n>\ 


_ 2(l-10n) 
' 3(l+u>) 


n > and n =£ -^ 


< n < jo 


n > To 


_ 2(l-12n) 
3(l+o>) 


n > and n / 75 


0<7J<i 


n > T2 



TABLE I: The relationship among r, oj and n, condition and candidate for late-time cosmic accelerated expansion in case 
/(G) = -^±L 



2G" 



From the above discussions, it is easy to sec that the results in the model are interesting. Compared with the f(R) 
models, f(G) models are even more complicated. Since the Hubble parameter can be expressed as H — j, the GB 
term turns into G = 



_ 24r J (r-l) 
t 



. If < r < 1, the early universe is in deceleration phase, which corresponds to the 

2 

3' 



matter dominated phase with r = ^, and if r > 1, the late universe is in acceleration phase 



Note that for the case of L m = p = top we can make similar discussions to ones in the case of L 
the same results as the ones shown in Table 1 due to the constant u). 



in 



-p, and obtain 



VI. CONCLUSIONS AND DISCUSSIONS 



In the present paper we have considered a modified f(G) gravity model with coupling between matter and 
geometry, described by the product of the Lagrange density of the matter and an arbitrary function of the Gauss- 
Bonnet term. The proposed action represents the general extension of the standard Hilbcrt action for the gravitational 
field, S = J d A Xy/— g{^ + [l + Xf (G)]L m } . The field equations and the equations of motion corresponding to this model 
show the non-conservation of the energy-momentum tensor, the presence of an extra-force acting on test particles and 
the non-geodesic motion. Moreover, in the modified f(G) gravity we have derived the energy conditions (SEC, NEC, 
DEC, WEC) when we consider L m — — p and L m = p, respectively. For the SEC and the NEC, the Raychaudhuri 
equation, which is the physical origin of them, has been used. From the derivation, we found equivalent results can 
be obtained by taking the transformations p — > p e $ ? and p — > p e $' into p + 3p > and p + p > 0. By means of these 
transformations, the DEC and WEC in the modified f(G) gravity have been also obtained. In order to exemplify how 
to use these energy conditions to constrain the modified f(G) gravity models, we have considered two specific models 
of f(G) gravity, i.e., fi(G) and /2(G) and given the corresponding constraints on the parameters in the fi(G) and 



12 

/2(G) models. Since there has been no reliable measurement for the snap parameter (s) up to now, we only focus 
on the WEC in this particular case. By analysis on Figs.l and 2 we have given the constraints on the parameters in 
the fi(G) and /2(G) models satisfying the weak energy conditions when L m = —p. By the similar discussions to the 
case of L m = —p, when L m = p the restrictions on the parameters n, b and a have been also illustrated in Figs. 3 
and 4, from which we have found that in the two different forms of L m the constraints on the parameters for the 
fi(G) model are nearly the same, but quite different for the /2(G) model. Furthermore, we have derived the stability 
criterion at the de Sitter point for the modified f(G) gravity models, which means that the modified f(G) gravity 
models may be stable. In addition, we have researched the conditions for late-time cosmic accelerated expansion in 
the modified f(G) gravity. Concretely, for the two different forms of L m , the relationship among r, u> and n have 
been respectively given in the model f(G) = — on^+b ' anc ^ ^ v using the conditions of power-law accelerated expansion 
and the equation of state of matter less than — ■=, the constraints on the parameter n have been obtained, which are 
exactly the same in the two different forms of L m . The candidate for late-time cosmic accelerated expansion would 
be either the effective quintessence (—1 < w < —h) or the effective phantom (w < —1), which could be determined by 
choosing n properly. Of course, other forms of f(G) gravity models with curvature-matter coupling will be considered 
in our following investigations. 
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